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ABSTRACT 

A new supersymmetric model for electrons with generalized hopping terms and 
Hubbard interaction on a one-dimensional lattice is solved by means of the Bethe 
Ansatz. We investigate the phase diagram of this model by studying the ground 
state and excitations of the model as a function of the interaction parameter, elec- 
tronic density and magnetization. Using arguments from conformal field theory we 
can study the critical exponents describing the asymptotic behaviour of correlation 
functions at long distances. 



PACS-numbers: 71.27. +a 75.10.Lp 05.70.Jk 



*e-mail: bed@itp.uni-hannover.de 
t e-mail: f rahm@itp . uni-hannover . de 



1 Introduction 



In recent years studies of one dimensional models of electronic systems have been a primary 
source to gain understanding of correlation effects in low dimensional systems. In particular 
the growing number of exactly soluble models such as the Bethe Ansatz integrable Hubbard 
and supersymmetric t—J models and their extensions have provided new insights into ground 
state properties of these systems Jl|] — ||. 

Different sources of interaction have been studied in these models: Apart from the influence 
of the on-site Coulomb repulsion (which is the main physical motivation leading to the Hubbard 
model) and the antiferromagnetic coupling of electrons leading to spin fluctuations (as present 
in the t—J model) the kinetic energy can been modified to include interaction effects. Such bond- 
charge repulsion terms reflecting the dependence of nearest neighbour hopping amplitudes on 
the occupation of sites affected were first discussed in Ref. ||. There have been extensive 
studies of the relevance of such additional interaction terms for example in their relation to the 
possibility of superconductivity based on electronic correlations (see e.g. (7], ||). Furthermore, 
several exact solutions for one-dimensional models of this type have been found (see e.g. f|, [|, |J). 

In this paper we consider a new integrable model containing generalized hopping integrals 
that has been found recently |ID, |TTJ. The Hamiltonian is given as 

Ti = - 53 ( c L c i+ia + h.c)j (t - X {rii- a + rii+x-a) + X m - a n i+ i - a ) 

- *3 J2( c l+i,A+u Ci i c ^ + h - c ) + u J2 n ^ n n ( L1 ) 

i i 

i A i 

In addition to the usual single particle hopping amplitude to and the on-site Coulomb integral 
U it contains the bond-charge interaction X, an additional coupling X correlating hopping 
amplitudes with the local occupation, and a pair hopping term with amplitude £3. In addition, 
the Hamiltonian contains coupling to a chemical potential \i and magnetic field h controlling 
particle density and magnetization of the system, respectively. 

For later convenience we introduce a different parametrization of the hopping integrals by 
X = t — t±, X = to — 2ti + ti- Studying the two particle scattering matrix S one finds two 
possible choices of the parameters tj and U where S satisfies a Yang Baxter equation resulting 
in candidates for models ( |1.1| ) that might be integrable by means of the Bethe Ansatz: first, 
choosing t = ti = t 2 (which implies X = X = 0) and t$ = the Hamiltonian reduces to the 
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well known Hubbard model [p]]. Another family of such models arises forQ 

-U = -t 3 = ±(to -h)^0, (ti) 2 = t t 2 (1.2) 



In an independent approach, the integrability of the model (|1 . 1| ) with ( |L2|) has been proven in 
the framework of the Quantum Inverse Scattering method where the Hamiltonian has been de- 
rived from a solution of the Quantum Yang-Baxter equation invariant under a four dimensional 
representation of gl(2\l) [|TIJ which is the symmetry underlying the (Bethe Ansatz soluble) 
super symmetric t-J model. 

Our paper is organized as follows: In the following section we shall discuss the symmetries 
of the Hamiltonian ( |1 . 1| ) at the integrable point (|1.2p . It turns out that there are two physically 
different regions to be studied corresponding to t > t 2 and t < t 2 (or positive and negative U), 
respectively. In Section 3 the Bethe Ansatz equations determining the spectrum of the model 
are derived. In Section 4 ground state properties and the spectrum of low-lying excitations at 
temperature T = are determined and in Section 5 we shall study finite size corrections of the 
spectrum to discuss the asymptotic behaviour of correlation functions. In the Appendix we 
discuss the completeness of the solutions obtained from these equations for small systems. 



2 Symmetries 



Owing to various symmetries of the Hamiltonian (LI) only the upper sign in the relation ( |1.2|) 
with positive t and t 2 has to be studied: 



To see this, we first note that the sign of ti is not fixed by the conditions ( |1.2|) . In fact, the 
unitary transformation 

C ia -> Cj CT (l - 2Tli - a ) (2.1) 

has the only effect of changing ti — > —t\. 

A particle-hole transformation performs a mapping between t ,t 2 > and t ,t 2 < (as a 
consequence of ( |L2p t an d t 2 necessarily have the same sign!): 

Tf. c la ^cl, (7=1,1 • (2.2) 

Applying this transform to the Hamiltonian we obtain (the irrelevant change of sign in ti is 
suppressed) 

H(t ,t 2 ,U = ±2(jt -t 2 ),fi, h) 

— > H(-t 2 , -t , U = ±2(t - h), n', -h) + (fjf — fx) L . (2.3) 



x The special cases t 2 =to/2 and t 2 — 2t have been discussed before in [JL2| 



Here // = 2 (to — t 2 ) — /■< and L is the number of lattice sites. 
The transformation 

T 2 : c l<T -> (-1)*^ , <r =t,| (2.4) 
changes the sign of the single particle dispersion resulting in 

H(t , t 2 , U = ±2(t - t 2 ), /x, h) -> H(-t , -t 2) U = ±2 (t - t 2 ), M, • (2-5) 
Applying both T\ and T 2 the sign in the first of Eqs. (|1.2j) is reversed (see Fig. Ill) 



H{t ,t 2 ,U = ±2{t -t 2 ),v, h) 

— > ft(t 2 , t , ^ = T2(t 2 - t ), /i', -/i) • (2.6) 

Hence t 2 and t are interchanged and at the same time the electronic density is changed from 
n e to 2 — n e . As will be seen later the Bethe Ansatz solution in the region to,t 2 > extends 
throughout the interval < n e < 2. Hence it is sufficient to consider the model with U = 
+2(t — t 2 ) in this region. 

As mentioned above the model can be constructed in the framework of the Quantum Inverse 
Scattering method based on a irreducible representation of the algebra g/(2|l). This is reflected 
in additional invariances of the Hamiltonian: apart from the SU(2) spin and number operator 

L L 
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=1 »=l 



^ e = £(n iT + n a ) (2.7) 
t=i 

which commute with the Hamiltonian for vanishing magnetic field there are four additional 
super symmetric generators fTt| , namely 



Qo = £(-1) V(r(! - ^,-) + ^,-) , ^ =t> 1 (2-8) 

i=l J o 

and their Hermitean conjugates Q\ satisfying commutation relations 

{Q T ,g i } = 0, Ql = 0, [H,Q a ] = (fi-2t + ah)Q a . (2.9) 



Fixing the potentials to = 2t , h = one obtains the supersymmetric model of Ref. |T0[ (up 
to the unitary transformation T{T 2 ). 

It is important to identify the full symmetry of the model since it is well known that the 
Bethe- Ansatz states are all highest weight states in this algebra and hence not complete flB] , i.e. 
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S + \^Bethe) = = Q a \ ^ Bethe) ■ Only after complementing the Bethe Ansatz states with those 
obtained by successive application of S ± and Q\ the complete set of eigenfunctions is found. 
We shall come back to this question at the Appendix. Note, that as a consequence of ( |2.9|) 
the number of particles in the states belonging to one gl(2\l) multiplet range from the number 
Nba m the Bethe-Ansatz state to Nba + 2. Hence, in the thermodynamic limit investigated 
below particle densities < n e < 2 can be studied directly. 



3 Bethe Ansatz solution in the thermodynamic limit 



Despite the derivation of the Hamiltonian ( |1 . 1|) in the framework of the Quantum Inverse 
Scattering method the spectrum of the model (which is obtainable in principle by means of 
the algebraic Bethe Ansatz) has not been found in |TIJ. The difficulty here is the complicated 
representation theory for the superalgebra gl(2\l). On the other hand, it is straightforward 
to determine the spectrum using the coordinate Bethe Ansatz: for models possessing internal 



symmetries as the one considered here the Schrodinger equation is solved with the Ansatz JL4 

*(X Q ) ^^A ai ,..., aN (P\Q)e^U^k Pj xA (3.1) 

where Q = {q±, . . . , q^} and P = {px, ■ ■ ■ ,Pn} are permutations of the integers {1, . . . , N} and 
Q is chosen such that Xq = {x qi < x q , 2 < . . . < x qN }. The coefficients A(P\Q) from regions 
different than Xq are connected with each other by elements of the two particle iS-matrix 

(P12 is a spin permutation operator). Here the charge rapidities dj are related to the single 
particle quasimomenta kj by ■#(&) = |tan(/c/2) and the dependence on the system parameters 
(|1 . 2|) is incorporated in the parameter c = (t Q — t 2 )/t 2 (varying in the intervals — 1 < c < and 
< c < oo). The A(P\Q) are determined in a second Bethe Ansatz for an inhomogeneous six 
vertex model resulting in the Bethe Ansatz equations (BAE) 

A9. _ i\ L m a. _ \ , -c 



= -n^f^. a=l,...,M (3.3) 

j=l A a — Vj — l 2 « =1 -A Q — A/3 — IC 



The length L of the system is assumed to be even and N e and M are the numbers of electrons 
and spin-J, electrons, respectively. Given a solution of ( |3.3j ) the eigenvalue of ( |1.1| ) in the 



corresponding state is 

/ AT \ N± 1 

E=(JMo-^ e -/i(-S-M)-toE-^XT- 

Solving (|3.3| ) one has to distinguish the two different regions c > and c < as the character 
of the $ — A solutions in these two cases is completely different. Note that the sign of the on-site 
Coulomb coupling U is the same as that of c because of ( |1.2| ) . Hence the situation is very similar 
to the Hubbard-model. Below we introduce some functions and their Fourier transforms that 
will be used in the following (y > 0): 



— f°° dk e-^e-m = — 11 



2tt 7-oc 2tt x 2 + y 2 /4 



1 f oc ^ 1 1 

Sy ^ ~ 2^J~oo dke 2cosh(?/fc/2) ~ 2y cosh(nx/y) ^ 3 ' 5 ^ 



i^x) = (a y *s y ){x) 



—Re\ib(l + i—) -ib(~ + i—)\ 



27n/ 

where (a * 6) (a;) = J dz a(x — z)b(z) denotes a convolution and r/> is the Digamma-function. 
3.1 Repulsive case (c > 0) 



In this case the solutions of (|3.3| ) consist of real i?j while the spin rapidities are known to be 
arranged in bound states of uniformly spaced sets of complex \ a , so called n- strings: 

\ n a 3 =\ n a + i(n + l-2 3 ) C - j = l,2,...,n (3.6) 

In the thermodynamic limit (L — > oo with particle density and magnetization ¥- being fixed) 
the solutions of the BAE (|3.3|) can be described in terms of densities p(i9) for charge rapidities 
and ph($) for the corresponding holes. Similarly, one introduces density distributions a n (aVi,fe) 
for the n-strings of spin rapidities (and corresponding holes). Using standard procedures one 
obtains the following system of coupled linear integral equations from the BAE ( |3.3| ) 

P + Ph = ai + Rc* P - s c * a 1;h 
cri+cri,h = s c *(a 2 ,h + p) (3.7) 
°n + °n,h = s c * (cr n +i,h + °n-l,h) , n>2 

The intervals in which the densities are nonvanishing depend on the state considered. The 
particle density and magnetization is related to p and a n through 

n e = -j- = I dd p{d) , 



1 (N e . . A n 
m * = L U 



oo 

oo 



-M) = ^-E / d\na n (X) 

' 1 n=l J -°° 



The energy density follows from (|3.4j) 

E 



e = — = (2t - A*K - hm z - 2nt / d$ ai(i?)p(t?) . (3.8) 

Li J — OO 

The equilibrium distribution functions p and a n have to be determined by minimization of 
the free-energy functional, T = E — TS, with the combinatorical entropy S of a particle and 



hole densities <5(A) and 8h(X) given by [15 

-L= d\ {(8 + S h )ln(5 + 8 h )-5\n(5)-8 h \n(5 h )} . (3.9) 

Li J — OO 

Introducing the following functions 

e e {0) = Tln(^) , e„(A) = Tln(^) (3.10) 

and considering p and a n ^h as independent functions we obtain by variation of T the following 
nonlinear integral equations for the functions e a 

e c = 2t — p — 2-rct ai + Ts c *\n(n(ei)) + TR c *\n(n(—e c )) 

£l = -Ts c *ln(n(e 2 )) +Ts c *ln{n{-e c )) (3.11) 

e n = -Ts c * (ln(n(e n+ i)) +ln(n(e n _i))) , n > 2 



with the distribution function 



n(e) = (l + e*) 1 . (3.12) 



Eqs. (|3.11|) have to be solved with the asymptotic boundary condition: 

\im— = h (3.13) 

n— »oo fi 

The free energy density is given by 

/•oo 

f = T d-d ai(i?)ln(n(-e c (i?))) . (3.14) 



This shows that the functions e a are to be identified as renormalized ("dressed") energies of 
the single particle excitations in the system. 

3.2 Attractive case (— 1 < c < 0) 

In this regime one has — in addition to the real charge rapidities and strings of spin rapidi- 
ties considered in the repulsive case — pairs of complex conjugated flf coupled to a real Xj as 
solutions of the BAE (O) 



%\c 



■df = A; ± (3.15) 
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Note that for c — > — 1 ( p.3|) would coincide with the BAE of the t-J-model obtained in ||, 
however this value is out of the range accessible for this model. Following the same programm 
as in the repulsive case we obtain in the thermodynamic limit 

P + Ph = S|c| * (?h + S\c\ * &l,h 

a' + a' h = ax_| c | + R\ c \ * o' h - s\ c \ * p 
ci + v\,h = S| C | * (a 2)h + p) 

<?n + <?n,h = S| C | * (a n+ i ih + a n - 1;h ) n>2 (3.16) 



where cr'(A) and cr' h (X) are the distribution function for the paired rapidities ( |3.15|) and corre- 
sponding holes. 



Particle density and magnetization of the state corresponding to a solution of ( |3.16D are 
given by the follwoing expressions: 



n e = d$ + 2 / dX cr'(A) , 

J — oo J—oo 

m z = - dd p{$) -J2 dX na n(X) , 

Z J — oo n=l J 00 

and the energy density is 

/oo 
-oo 

-2nt d\ (oi +k | + Oi_ k |) (AK(A) . (3.17) 



Minimizing the free energy (with a' h as additional independent function) and defining the 
dressed energy of the paired rapidities as 



e p (X) = T\n[^) (3.18) 
we obtain the thermodynamic Bethe Ansatz equations for the attractive case 
e c = Ts\ c \ *ln(n(£i)) - Ts\ c \ * ln(n(e p )) 

e p = 4t - 2/i - 2irt (a 1+ \ c \ + ai_| c |) - T (a 2 | c | * ln(n(-e p )) + a\ c \ * ln(n(-e c ))J 
ei = -Ts| c| *(ln(n(e 2 ))+ln(7i(-e c ))) 

e n = -Ts| c |*(ln(n(£ n+ i))+ln(n( £n _i))) , n>2 (3.19) 
to be solved with the field boundary condition ( |3.13| ). The free energy density is given by 

/ = T dtf (ai(tf) ln(n(-£ c )) + (ai_| c , + a 1+ , c ,) ln(n(-£ p ))) . (3.20) 

t/ — oo v ' 



4 Ground state and excitations at T = 



Now we want to examine properties of the zero temperature ground state for the two cases. In 
this limit the distribution function n ( |3.12| ) in the thermodynamic BAE reduces to 

limTln(n(5)) = -5 + (4.1) 

where 5 + > and 5~ < are the positive and negative parts of the function 5 = 5 + + 5~, 
respectively. At the same time it is clear that the ground state configuration corresponds to 
the filling of all states with negative dressed energy e a . 



4.1 Repulsive case (c > 0) 



From (|3.11|) we find that £ n >i(A) > for all A. Using the asymptotic condition ( |3.13 ) we obtain 
from (1DID with 



^ 2t - p - 27rt ai - | 



V 



h 



V 


( a c \ 


■ 


{■'■-) 











(4.2) 



As in (T^l one can prove that e c (^) and S\{X) are monotonically increasing functions of \&\ and 
|A|. Consequently, they are negative in the intervals [— Q,Q] and [— B,B]. For h = two 
possible ground state configurations are to be considered: the ferromagnetic state (M = for 
n e < 1 and M = N e — L for n e > 1) and the antiferromagnetic one (M = N e /2). The energy 
of the ferromagnetic state at fixed density n e is simply 



- sin irn e 

IT ° 



for n P < 1 



< r.vi = -2/,, { (4.3) 

sin7rn e - c(n e - 1)) for n e > 1 

For the antiferromagnetic state one obtains from ( |3.7| ) that it corresponds to a filled band of 
1-strings, i.e. B = oo. Hence <j\ can be eliminated by Fourier transform and ( |3.7|) simplifies to: 

'■Q 



d-d' R c {d - 0')p(&) 



(4.4) 



Varying Q one obtains any filling between n e = and n e = 2: For small Q ( |4.4|) can be solved 
by iteration and for Q — > oo using Wiener Hopf techniques |l7j with the result 



n,, 



f + 8^g 2 + 0(Q 3 ) , forQ^O 
2-^(1 + ^)+°^) forg^oo. 



(4.5) 
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In the low density limit we find that the ground state of the system is indeed antiferromagnetic. 
The energy difference to (|4.3|) is 



^nt + Oinl) , for n e -> 



e F M ~e = { (4.6) 
^(2 -Tie) +o((2-n e ) 2 ) forn e ^2. 

In Figure ^| we present numerical data for the dependence of the (antiferromagnetic) ground 
state energy of the system as compared to the ferromagnetic one for various values of the 
parameter c.f\ 

In the free fermion limit c — > (|4.4|) simplifies to 



p (f) = a 1 (&) + - I' d&' <5(tf - tf')p(tf') ( 4 -7) 



2 V- . 

and the ground state energy is the expected result for this system 

e = sm(— ) . (4.8) 

7T Z 

In the strong coupling limit c — > oo corresponding to ti = t 2 = 0, i 3 = — 1 , U = 2t , the 
groundstate is degenerate with the ferromagnetic state (|4.3| ). 

There are two types of excitations that are to be considered in the low energy sector: first, 
there are objects carrying charge ('holons') corresponding to particle or hole like excitations 
in the ground state configuration of charge rapidities with energy |e(i?)|. Furthermore, there 
are spin carrying objects ('spinons') corresponding to holes in the distribution of real spin 



rapidities. From ( [4.2|) their energy is found to be 

u ; 2cJ-q coshf(A-^) v ; 

The physical excitations (for even particle number N e ) are even numbers of these objects 
forming a continuum of spin waves without gap. The energies of the spin rapidity strings of 
length n > 1 vanish. 

Increasing the magnetic field the magnetization grows until it reaches its saturation value 
| at the critical field h c . For h = h c the interval for the A-integration vanishes, i.e. B = 
(corresponding to £i(A = 0) =0). For < n e < 1 we find 

8*o c(4Q 2 + 1) arctan(2Q) - (4Q 2 + c 2 ) arctan( 22) 

hc = T FipHlj (4 - 10) 

with Q = |tan( 2E | £ ). In the limiting cases considered above this expression becomes 

7TTI 

lim h c = , lim h c = 4t 2 , \imh c = 4t sin 2 ( — -). (4.11) 

Q— >0 Q^oo c^O 2 



2 In Rcf. [12| the ground state is claimed to be /erromagnetic in this regime for densities n e < 1. As is clear 



from Eq. (4.6), this is not correct. 
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4.2 Attractive case (— 1 < c < 0) 



Due to ( p.!9|) the dressed energies of the spin rapidities £ Tl >i(A) are always positive. Performing 



the limit T -»• in j ^Wjj with Eqs. Q and ( P7TB|) we obtain 



V 



^ «2|c| «|c| 



£ P 



y 0| c | 



(4.12) 



4t - 2/x - 27rt («i+|c| + ai-|c|> 
y 2t — fi — 2-nt a ai — | 

As in the repulsive regime one can prove that e c {d) and e p (X) are monotonically increasing 
functions of the modulus of their arguments. Hence they are negative in the regions [-Q, Q] 
and [-B, B\. Again we find that the ground state of the system is antiferromagnetic for h = 
(see Fig. |3[). In this regime the ground state configuration consists of paired rapidities only 
(Q = 0). Their density is obtained from ( |3.16| ) which simplifies to 

or'(A) = a x _| c |(A) + Oi+| c |(A) - / dfi a 2 | c |(A - n)a (fi) (4.13) 

J — B 

Again, this is the ground state configuration for any filling < n e < 2 since 

ff T^ + 0(5 2 ) for B^O, 

n e = (4.14) 

[ 2 -^(l-J^)+0(i) forB-oo. 
For c — > (|4.13|) simplifies to: 

or'(A) = 2ai(A) - [ B d[i 5(X - ^)a'U) (4.15) 

J-B 

and we obtain for the ground state energy ( |4.8| ) of the free fermion system. 

From Eq. ( |4.12j ) we obtain the dressed energy for excitations corresponding to real charge 
rapidities: 



e c (#) = -- + _/_ rfA _^„_ Q 7^373 e p (A) (4.16) 



/i 1 

2 + ]c\Jb "" cosh(($ - A)tt/c) 
Note that for vanishing magnetic field there is a gap A c (n e ) = £ c (0) for the creation of unpaired 

electrons. A c is a monotonically falling function of n e with its maximum at A c (0) = > 
(see Figure f|). The only massless excitations in this regime are charge density waves corre- 
sponding to excitations within the band e p (X). 

In an external magnetic field the nature of the excitations in the system change: For h c \ = 
2A C the gap for charge excitations closes and for h C 2 the system undergoes a transition into 
the saturated ferromagnetic ground state. The latter corresponds to B = (or, equivalently 
£T P (0) = 0). For < n e < 1 we find: 

8 t oC 2 8t (4B 2 + c 2 )arctan(^)-|c|(4 J B 2 + l)arctan(2 J B) 
kc2 -\^ + ~ (c 2 -l) (41^ + 1) (4 ' 17) 
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with B = |tan( IL ^-). We obtain the following limiting cases: 



lim h c2 = At 2 lim h C2 = 2A C (0), \imh C2 = 4t sin 2 (— -). (4.18) 

B^oo B— >0 c^O 2 

The nonvanishing of h c2 in the low density limit is a direct consequence of the gap of e c . 



5 Finite size corrections and critical exponents 

We now want to study the finite size corrections of the spectrum to discuss the asymptotic 
behavior of correlation functions. Again the repulsive and the attractive case are completely 
different and have to be treated separately. 

5.1 Repulsive case (c > 0) 

As found in Section 4.1 above for T = the ground state and low lying excitations are obtained 
from solutions of the BAE (|3.3|) with real i?'s and A's. Hence, we have the same situation as in 
the repulsive Hubbard model and following the procedure in Ref. we obtain the finite size 
corrections of the ground state energy as 

E -Le = -^(v c + v s ) + o(-\ , (5.1) 

where v c and v s are the Fermi velocities of charge and spin density waves, respectively: 

1 <(Q), ".= n ■ (5-2) 



2np(Q) cv ^ y ' 2na l {B) 
Similarly the energies and momenta of the low lying excitations are given by 



2ir r 

E(AN, D) — Le = — {v c (A+ + A;) + v s (Aj + A s 



+ o 



2lX r i 

P(AN,D)-P = -j-[At-A: + At-Aj}+2D c P F ^ + 2(D c + D s )V F , l (5.3) 

with the conformal dimensions 

+ / n ( Z SS AN C - Z CS AN S \ 2 
2A±(AN,D) = iz cc D c + Z sc D s ± ss ^7^y -J + 2N t 



c 

2 



2A±(AN,D) = [Z CS D C + Z SS D S ± ^f^)^ ) + ^ ' (5 ' 4) 

and the Fermi momenta Pf,t(!) = f K±2m z ) of spin up (down) electrons. The elements of the 
two component vectors AN and D characterize the excited state: AN has integer components 
denoting the change of the number of electrons and down spins with respect to the ground 
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state. D c and D s describe the deviations from the symmetric ground state distributions. They 
are integers or half-odd integers depending on the parities of AN C and AN S : 



AN C + AN S 



mod 1 , -Do 



AN r 



mod 1 



(5.5) 



(5.6) 



The matrix 

UQ) U(Q) 

Z sc Z ss J \ £cs(B) £ss(B) 

parametrizing the conformal dimensions ( |5.4| ) is given in terms of the so called dressed charge 
matrix which satisfies a linear integral equation similar to (|4.2|) for the dressed energies 




(5.7) 



As shown above the ground state at vanishing magnetic field corresponds to B = oo and 
with the aid of the Wiener-Hopf method [pl| (|5.6|) simplifies to 

^ /,.... /,... I I C ..((./I u I 

(5A 



V 



UQ) o 
\UQ) 4 



2 scv^y 2 

where £ c is defined as the solution of the following scalar integral equation 

rQ 



(5.9) 



As for the density one can solve fl5.9|) near Q = and Q = oo with the result 

1 + ^Q + 0(Q 2 ), forQ^O 
V^(l-^)+ (i) forQ^oo. 
Hence the range of variation for the exponents determining the long distance asymptotics of 



UQ) 



(5.10) 



the equal time correlators is the same as in the Hubbard- and t—J- model p9| , p20| . Introducing 
9 = 2^(Q) the singularity of the momentum distribution function at the Fermi point is found 
to be 



n a (k) ~ J dx e- ikx (c x ^(t = + )c , a (t = 0)) 
oc sgn(k - V F )\k - V F \ U , 



v 



1 9 1 
9 + 16 ~ 2 



(5.11) 



with a variation of the exponent v in the interval < v < | which shows the expected Luttinger 
liquid behaviour of this system. 
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Similarly, we obtain for the density density and singlet pair correlation functions (Vf,i = 



G nn (x) = ( {n x] + 71,; ) (n 0T + n |)> 

~ n 2 e + A 1 cos(2V F x + ip 1 )x' {1+9/4) + A 2 cos(AV F x + (p 2 )x- d + A 3 x~ 2 (5.12) 
G®(x) = (ci +liTC i aCuCo , T )~Acos(2P F a ; + ^x-( 4 /^/ 4 ) 

The leading order of the density-density correlator is given by the A\ term with 3/2 < 1 + 0/4 < 
2. Comparing this with the leading term of the singlet-pair correlator 5/2>4/0 + 0/4>2we 
see that density fluctuations are dominant. 

In Figure [| we show lines of constant £ C {Q) (hence identical critical behavior) in the n e -c 
parameter plane. Note that the strong coupling result £ C {Q) = 1 is found for less than half 
filling only. Beyond half filling the density dependence of the dressed charge is for c = oo 

n e for n e — > 1 

ic={ _ (5.13) 
y/2(\-\(2-n e )) forn e ^2. 



As in |21| f° r the Hubbard model this analysis of the critical behaviour can be extended to 
the case of magnetic fields. For small fields h < h c one has to expect logarithmic singularities 
in the exponents while for fields h > h c the ground state is a saturated ferromagnetic one and 
spin density waves become massive giving a scalar dressed charge instead of ( |5.8| ). 

5.2 Attractive case c > 

As discussed in Section 4.2 for T = and h < h c i there is only one branch of massless excitations 
within the band e p ^\ The finite size corrections to the energies of the low lying excitations are 
given by 

E(AN P ,D P )-Le = + A p ) + o Q 



2-7T 

P(AN P ,D P )-P = T (A;-A p ) + 2D P V F (5.14) 



with 



2A±(AiV p , D p ) = (z p (B)D p ± + 2N+ (5.15) 



3 In the analysis of the asymptotics of correlation functions for the model with c = — 1/2 in j23] the existence 
of a second branch of massless excitations in the band of real charge rapidities e c is assumed. However, as 
shown in Sect. 4.1 these have a gap for h < h c \. Hence the results in [E2| are incorrect. 
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and charge density wave velocity v p = e' (Q)/(27ra'(Q)) The dressed charge £ p is given by 

&(A) = 1 - f B d\'a 2lcl (X - A')£ P (A'). (5.16) 

J — B 

With the same techniques as above we obtain 

1-^B + 0(B 2 ), forfi^O 
UB) = { (5.17) 

The leading terms in the asymptotics of the equal time correlators as a function of 9 = 2^(B) 
are the same as in the (attractive) Hubbard model |23| 

G nn {x) ~ nl + a ^W**) + ^ , G (o) (x) „ x -i/* (518) 

Comparing the leading exponents of these two correlators we see that the correlation of pairs 
(1/2 < 1/8 < 1) overwhelms the density-density correlator (2 > 9 > 1) fo arbitray n e . So as 
in the attractive Hubbard model [p3| we can conclude that the particles are confined in pairs 
which is reflected in the structure of the Bethe Ansatz ground state configuration. 

In Figure ^|we show lines in the n e -|c| parameter plane with identical critical behavior. 

For h > h c i charge and spin excitations are massless and the dressed charge is a 2 x 2 matrix 
as in the repulsive case. The same situation occurs in the attractive Hubbard model 
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A Completeness of the Bethe Ansatz states 

As mentioned above the Bethe Ansatz states do not form the complete set of eigenstates of the 
systen ( |1.1|) but are the highest weight states of the g/(2|l) superalgebra. Complementing the 
Bethe Ansatz states with those obtained by the action of the gl(2\l) shift operators one obtains 
additional eigenstates. The completeness of this extended Bethe Ansatz has been proven (based 
on a string hypothesis (|3.6|) for the solutions of the BAE) for some models such as the the spin 
I Heisenberg chain, the supersymmetric t—J model and the Hubbard model J13, |24|, [25|]. In this 



appendix we present the study of the completeness for the two-site system together with some 
remarks on L > 4. 
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New eigenstates of the system are generated from the Bethe Ansatz states by acting with 
the total spin operators S~ , S + and the supersymmetry generators Q a , Q\. As a consequence 
of the anticommutativity of the latter (|2.9| ) the resulting multiplet contains states in the N e -, 
(N e + 1)- and (N e + 2)-particle sectors which are (we suppress the spin- multiplicity): 

l W < ( ) Sk |* e3) . (a.i) 

[ *Q2> J 

As shown above, the ground state of the model for fixed number of particles is always a spin 
singlet. As a consequence of ( |A.1|) it is member of a gl(2\l) quartet, the same situation as in 
the related supersymmetric t—J model 



IV*> = 


\N e 


= 0, 


M 


M = 


\N e 


= 1, 


M 


M = 


\N e 


= 2, 


M 



Solving the BAE (|3.3|) in the simplest case of the L = 2 system we obtain three regular 
Bethe Ansatz states \ipi) with energy Ei (at the supersymmetric point \i = 2t Q , h = 0): 

0) = |0) , Eh = Q 

0) = \k = 0, T) , E 2 = -4t (A.2) 

1) oc |V n > - |V; T > + rd^o) + l^02» , = -4(t + *a) 
with 

\^a 2 ) = 4,^4^)0) , I^q) = 4, T 4,i|0> » l^oa) = 4, t 4,J°> ■ ( A - 3 ) 
Regular Bethe Ansatz states are those corresonding to solutions of ( |3.3|) with finite d and A 

The one- and two particle descendants of are found to be the momentum 7r spin-doublet 
\k = n, a) and the spin singlet 

IM-M-T^ + IM ■ (A.4) 

Analogously we find the descendants of l^) to be the following (degenerate) triplet and singlet 
states in the two particle sector 

IV>n>> I'M - I'M ( A - 5 ) 

and the doublet of zero- momentum single hole states \kh = 0, a). Finally, \ips) leads to doublet 
of momentum ir hole states \kh = 7r, cr) and the completely filled state | "022 ) - 

Hence the Bethe Ansatz extended by means of the supersymmetry does indeed give the 
complete spectrum of states on the two site lattice. Note, that l^) is always the ground state 
of the two particle sector for the range of parameters considered here. The difference between 
the repulsive and attractive regime is the larger amplitude of the states containing local pairs 
in the latter. 
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For general L regular Bethe Ansatz states will exist for particle numbers up to 2{L — 1). 
Considering L = 4 as an example one has to find 35 regular solutions of the BAE to generate 
a complete set of eigenstates. Four of these are states with N e > L. 
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Figures 




Figure 1: 

Range of parameters to, t 2 for which the model ( |1 . 1| ) is integrable as a consequence of ( |1.2| ) for 
U = +2(t — t 2 ) (left) and U = — 2(t — t 2 ) (right). The dots mark the model introduced in 
Ref. 0. 
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Figure 2: 

Energy of the antiferromagnetic ground state of the system ( |1 . 1|) vs. electron density in the 
repulsive regime for various values of the reduced coupling constant c. For comparison, the 
energy of the ferromagnetic state (|4.3|) is also included. Note, that for c — ■> oo the ferro- and 
antiferromagnetic states are degenerate. 
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c=-0.1 



c=-0.5 



c=-0.8 



Figure 3: 

Energy of the antiferromagnetic ground state of the system ( |1 . 1|) vs. electron density in the 
attractive regime for various values of the reduced coupling constant c. 
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oo 




Figure 5: 

Contours of constant £ C (Q) (and hence identical critical exponents) in the n e -c parameter plane 
of the repulsive model. £ C (Q) varies between 1 (at low densities) and y/2 (the free fermionic 
case) for finite c. 
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|c| 




Figure 6: 

Contours of constant £ P (Q) (and hence identical critical exponents) in the n e -|c| parameter 
plane of the attractive model at small magnetic fields h < h c \. £ P (Q) varies between 1 and 
1 / \[2 for any finite c. 



23 



